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One of the great challenges in the ecology of infectious diseases is to understand what drives the emer¬ 
gence of new pathogens including the relationship between viruses and their hosts. In the case of the 
emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), several studies have shown 
coronavirus diversity in bats as well as the existence of SARS-CoV infection in apparently healthy bats, 
suggesting that bats may be a crucial host in the genesis of this disease. To elucidate the biogeographic 
origin of SARS-CoV and investigate the role that bats played in its emergence, we amplified coronavirus 
sequences from bat species captured throughout Thailand and assessed the phylogenetic relationships to 
each other and to other published coronavirus sequences. To this end, RdRp sequence of Coronavirinae 
was targeted by RT-PCR in non-invasive samples from bats collected in Thailand. Two new coronaviruses 
were detected in two bat species: one Betacoronavirus in Hipposideros larvatus and one Alphacoronavirus 
in Hipposideros armiger. Interestingly, these viruses from South-East Asia are related to those previously 
detected in Africa (Betacoronavirus- b) or in Europe (Alphacoronavirus & Betacoronavirus- b). These findings 
illuminate the origin and the evolutionary history of the SARS-CoV group found in bats by pushing for¬ 
ward the hypothesis of a Betacoronavirus spill-over from Hipposideridae to Rhinolophidae and then from 
Rhinolophidae to civets and Human. All reported Betacoronaviruses-b (SARS-CoV group) of Hipposider¬ 
idae and Rhinolophidae respectively cluster in two groups despite their broad geographic distribution 
and the sympatry of their hosts, which is in favor of an ancient and genetically independent evolution 
of Betacoronavirus-b clusters in these families. Moreover, despite its probable pathogenicity, we found 
that a Betacoronavirus-b can persistently infect a medium-sized hipposiderid bat colony. These findings 
illustrate the importance of the host phylogeny and the host/pathogen ecological interactions in the 
description and the understanding of pathogen emergence. The host's phylogeny, biogeography and 
behaviour, combined with already described roles of pathogen plasticity and anthropic changes are likely 
to be co-factors of disease emergence. Elucidating the common ancestor of Hipposideridae and Rhinolo¬ 
phidae is key to understanding the evolutionary history of actual betacoronaviruses and therefore to get 
an insight of the deep origin of SARS-CoV. 

© 2011 Elsevier B.V. All rights reserved. 


1. Introduction 

Pathogens mostly emerge in densely human populated areas as 
did the Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), 
which had the potential to cause the first viral pandemic of the 
XXIst century (WHO, 2003; Jones et al„ 2008). SARS-CoV most 
likely originated from bats through an intermediate species that 
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could be the masked palm civet Paguma larvata (Viverridae) or an¬ 
other bridge-host, itself very probably contaminated by a SARS-like 
CoV from horseshoe bats (Song et al., 2005; Li et ah, 2006; Cui et ah, 
2007). Several elements indicate that this bridge-host could pre¬ 
sumably be a Carnivora as suggested by both detections in wild 
animals from markets or experimental infections (Guan et ah, 
2003; Van den Brand et ah, 2008). Several other mammalian or 
avian species are susceptible to coronaviruses infection (displaying 
enteric, hepatic, neurological, or respiratory tropism), but little is 
known about their viral ecology and host interactions. Coronavi¬ 
ruses are members of the order Nidovirales which contains three 
families of single stranded positive RNA viruses: Arteriviridae, 














M. Ar Couilh et al. /Infection, Genetics and Evolution 11 (2011) 1690-1702 


1691 


Coronaviridae and Roniviridae. Coronaviridae are divided in two 
subfamilies: Coronavirinae and Torovirinae. Because of serological 
and genotypic differences, coronaviruses were traditionally 


divided into three groups labeled group 1, II and III (Gonzalez 
et al., 2003; Gorbalenya et al., 2004). These three groups have re¬ 
cently been elevated at the genus level and named respectively 


1 A 


AJ311317Avian IBV Beaudette CK 2001 
- EU095850Turkey CoV MG10 2008 


- EF065509 Bal CoV HKU5 1 China 2006 
.TEF065505 Bat CoV HKU4 1 China 2006 

• 00648 794 Bal CoV 133 China 2005 1 

- EF065513 Bal CoV HKU91 China 2006 
“ AB539081 Bat CoV Philippine Diim 2008 

• DQ415914 Hum CoV HKU1 N18genoA2006 

• FJ647225 Murine Hepatitis Virus A59 2009 



'0022305 Bat SARS CoV HKU3 1 China 2005 
, 'Q412043 Bat SARS CoV Rml China 2005 
^00412042 Bat SARS CoV Rfl China 2005 




:hma 2004 



Gammacoronavirus 


Betacoronavirus 


—Avian coronaviruses 


— Domestic animals coronaviruses 
—Chiropteran coronaviruses 

— Human coronaviruses 

—SARS coronaviruses 


0.0060 



DQ811789 TGEV virulent Purdue 2006 
DQ010921 Feline CoVFIPV7911461979 
AF516906 CanineCoV UWSMN 1 2002 

EF203065 Bat CoV HKU2HK46 China 2006 
GUI90238 Bat CoV BB98 41BGR Bulgana 20CM 
AF304460 Human CoV 229E 2000 


■ FJ710045 Bat CoV Hip Ghana Kwarr 

• ", ' 


ated 2004 


AY518894 Human Cov pneumoi 

AF353511 Porcine EDVCV777 1993 

EF544563.1 Rocky Mountain Bat CoV 11 US 2006 

E DQ648858 BtCoV5l2 China 2005 
AB539080 Bat CoV Philippine Dilim 2008 
DQ249224 Bat CoV HKU 6 1 China 2005 
EU375864 Bat CoV Pnat Germany 2007 
EU375866 Bat CoV M dau Germany 2007 
EU769557 Bat CoV Trinidad 2007 
108 ■ ■" - - - ,v ■' - Tf nil nr it '<§> 

j,GU190234 BatCoV BM48 39BGR Bulgaria 2008 
1 GU19C235 Bat CoV BM96 05BGR Bulgaria 2006 
EU834951 Bat CoV M. mac Australia 2008 
GUI 90239 Bat CoV BNM98 30BGR Bulgana 2008 
GUI 90240 Bat CoV BR98 55BGR Bulgana 2008 
EU420137 Bat CoV IB AFCD307 China 2008 
EU420138 Bat CoV 1AAFCD62 China 2008 
Ti DQ249226 Bat CoV HKU7 1 China 2005 
v GUI90243 BatCoV BR98 30BGR Bulgaria 2008 
:U834956 Bat CoV M. aus Australia1996 
EU834955 Bat CoV M, sch Auslraha 2007 
■ EU834953 Bat CoV R. meg Australia 2007 
EU420139 Bat CoV HKU8 AFCD77 China 2008 
iGU 190242 Bat CoV BR9818BGR Bulgaria 2008 
GUI 90246 Bat CoV BR98 40BGR Bulgana 2008 


Alphacoronavirus 


1 B 



Fig. 1 . Phylogenetic reconstruction of main coronaviruses based on the analysis ofof 333 bases of the RNA-dependent RNA polymerase protein (partial nspl2) and 
geographical localization of Hipposideros larvatus Betacoronavirus Tata and Hipposideros larvatus Alphacoronavirus Ratcha. (A) Bayesian phylogeny of coronaviruses showing 
host taxonomy of Betacoronavirus- b group and reporting two new coronaviruses (posterior probabilities are shown for each node). (B) Repartition of hosts families (Hutson 
2003), sampling sites and HI BCoV T & Ha BCoV R spatial position. 
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Alphacoronavirus, Betacoronavirus and Gammacoronavirus by the 
coronavirus Study Group of the International Committee for 
Taxonomy of Viruses, thus invalidating the "Coronavirus" genus 
name (http://talk.ictvonIine.Org/media/p/1230.aspx). As other Nid- 
ovirales, Coronaviridae such as SARS-CoV possess a mutation prone 
RNA-dependent RNA-polymerase (RdRp) associated with high gen¬ 
ome mutability (Castro et al., 2005). Furthermore, they frequently 
recombine in vivo (Lai et al., 1985; Makino et al., 1986; Keck et al., 
1987; Jia et al., 1995; Herrewegh et al., 1998; Lee and Jackwood, 
2000; Stadler et al., 2003). This genetic variability confers to coro- 
naviruses a high potential of evolution that allows them to some¬ 
times overcome species barriers and host specificity (Baric et al., 
1995, 1997) (Thackray and Holmes, 2004). From a more general 
viewpoint, these characteristics are known to account for the fre¬ 
quent involvement of RNA viruses in disease emergence (Domingo 
and Holland, 1997; Chomel, 1998; Daszak et al., 2000; Cleaveland 
et al„ 2001; Woolhouse and Gowtage-Sequeria, 2005). 

Until now, Asian SARS-related coronavirus (SARSr-CoV, from 
Betacoronavirus -b group - Fig. 1A) RNA has been detected in only 
one bat family: the Rhinolophidae Gray, 1825 (Table 2). Therefore, 
this monogeneric bat family is considered as the host at the origin 
of the deadly SARS-CoV outbreak (Lau et al., 2005). However, other 
bat families are involved in the maintenance and the diversifica¬ 
tion of a wide range of coronavirus species belonging to both 
Alphacoronavirus and Betacoronavirus genera (i.e. Vespertillionidae, 
Phyllostomidae, Pteropodidae, Molossidae, and Hipposideridae 
families have been positively tested for coronavirus nucleic acid) 
(Chu et al., 2006; Tang et al., 2006; Dominguez et al., 2007; 
Carrington, 2008; Brandao et al., 2008; Tong et al., 2009; Pfefferle 
et al., 2009). Whether or not another bat family could be involved 
(in addition to Rhinolophidae) in the natural history of betacornav- 
iruses would have been an interesting question but since the dis¬ 
covery of SARSr-CoV in wildlife, most of the attention and 
sampling efforts have been focusing on Rhinolophidae (Table 2). 
Several species from this family were shown to carry different 
SARSr-CoV such as SARSr Rhinolophidae Bat Coronavirus strain Rp3 
in Rhinolophus pearsoni and SARSr Rhinolophidae Bat Coronavirus 
strain Rfl in Rhinolophus ferrumequinum in the Chinese provinces 
of Guangxi and Hubei, respectively (Lau et al., 2005; Li et al., 
2005; Tang et al., 2006). The sampling pressure on both Rhinolo¬ 
phidae (1148 samples in 12 species - Table 2) and Vespertilionidae 
(900 samples in 17 species - Table 2) has induced an over repre¬ 
sentation of these families compared to others (Lau et al., 2005; 
Li et al„ 2005; Tang et al., 2006). By contrast, Hipposideridae 
(206 samples in 6 species) were under-represented in previous 
studies (Poon et al„ 2005; Woo et al., 2006). The higher northern 
distribution limit of Rhinolophidae (versus Hipposideridae) may 
partially explain the sampling bias in non-Southern Chinese prov¬ 
inces (Fig. IB). This bias illustrates the need for testing other sym- 
patric species for betacoronaviruses to study their natural history 
and emergence process. 

The biogeography and the natural history of the Betacoronavi- 
rus -b species, among which are SARSr-CoVs, remain poorly de¬ 
scribed. However, Betacoronavirus have been detected in China, 
Africa, Europe and even Australia (Guan et al., 2003; Lau et al., 
2005; Li et al., 2005; Muller et al., 2007; Tong et al., 2009; Quan 
et al., 2010; Drexler et al„ 2010) (Fig. IB). The biogeographical ori¬ 
gin of the ancestors of SARSr-CoVs (from Betacoronavirus -b sub¬ 
group) is therefore far more extended than the epicentre of the 
human epidemic (Southern China). Understanding the biogeogra¬ 
phy of this group would be essential for assessing the risk for 
emergence of another SARS-CoV like agent. South-China and its 
Pearl River delta combined various ecological and anthropic factors 
that have led the SARS-CoV to emerge but the distribution of the 
viral clade at the origin of SARS-CoV is much wider than this region 
and many locations still remain to be investigated. 


Emergence results from a combination of numerous factors. 
Environmental changes (ecological, climatic and agricultural), 
whether they are man-driven or not (Schrag and Wiener, 1995; 
Kuiken et al„ 2003), globalization (trade or travel) and societal 
changes (land-use, demography), are among the most reported 
ones (Reperant, 2010), but others have been highlighted recently 
like the hosts’ phylogenetic relationships (Daszak, 2010; Streicker 
et al„ 2010), the variation of biodiversity (Keesing et al., 2010) 
and species interactions in ecological communities whether be¬ 
tween hosts or between hosts and pathogens (Cleaveland et al., 
2001; Woolhouse and Gowtage-Sequeria, 2005; Reperant, 2010). 
Consequently, the biogeography and host-pathogen relationships 
of Betacoronavirus -b are pivotal to understanding their radiation. 
The diversity of the SARSr-CoV has been described by an increasing 
number of published studies and the number of discovered and se¬ 
quenced viral taxa has improved the accuracy of the phylogenetic 
reconstruction of this group. 

In this study we explored the biodiversity and the ecology of 
Coronavirinae, enriched the current phylogeny and challenged the 
biogeography of this clade. More particularly, we focused on the 
ecology of the Betacoronavirus -b sub-group (Fig. 1A) by testing 
hypotheses on the origination of these viruses not only in Rhinol¬ 
ophidae but also in their sister taxa Hipposideridae. Considering 
that Hipposideridae might host a diverse community of coronavi- 
ruses comparable to that of Rhinolophidae because of their phylo¬ 
genetic proximity and comparable diversity (Nunn et al„ 2004; 
Simmons, 2007), it is very probable that the prevalence and diver¬ 
sity of coronaviruses hosted by this family has been underesti¬ 
mated because of biased sampling (Table 2 - Lau et al., 2005; Li 
et al., 2005; Tang et al., 2006). We hypothesize also that some 
bat taxa could play a refuge role for betacoronaviruses and there¬ 
fore constitute a source of infection for other species, consistently 
with the recent application of the theory of island biogeography 
(MacArthur and Wilson, 1967) to emerging pathogens (Reperant, 
2010). This study is particularly focused on, but not restricted to, 
Betacoronavirus in Hipposideridae in South-East Asia (Table 1 ). 

2. Materials and methods 

To test our hypotheses, we targeted, by RT-PCR, the Coronaviri¬ 
nae RdRp RNA sequence (more precisely a coding fragment of a 
subunit of the RdRp, the non-structural protein 12 - nspl2 -) in 
bat samples and compare it with the published sequences down¬ 
loaded from public online databases. The SRAS-CoV nspl2 and 
nspl3 exhibit the highest conservation profiles among proteins 
from all coronaviruses (Stadler et al., 2003). Moreover, as RdRps 
(nspl2) are involved in replication and transcription processes 
and therefore play a central role in the virion offspring genome 
synthesis (the vertical information transmission), its relatively 
conserved ribonucleic acid sequence is appropriate for both viral 
detection and phylogeny reconstruction. 

2.1. Sampling 

From 2006-2008, samples (feces, urine and oral swabs) were 
collected from bats at 26 sites in Thailand (Fig. IB and Table 1). 
Sampling locations were selected considering land cover and geo¬ 
logical data as forests are known to enhance biodiversity (Harvey 
et al., 2006; Trisurat et al., 2010) and limestone hills harbour 
numerous caves. The rocky forest-covered hills surrounded by cul¬ 
tivated lowlands, considered as potential refuge areas, were partic¬ 
ularly sought. Sampling sessions covered an area of approximately 
800 km from North to South Thailand by 300 km from West to 
East. Dates and locations are described in Table 1 and Fig. IB. Sites 
were sometimes sampled several times to confirm the results, to 
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Table 1 

Detection of coronaviruses in saliva (Oral Swab), feces and urine in bats from Thailand. 


Date 

Family 

Genus 

Species 

O Swab 

Faeces 

Urine 

Province 

Cave 


N 

+ 

N 

+<%) 

N 

+ 


2008-01 

Craseonycteridae 

Craseonycteris 

thonglongyai 

- 

- 

3 

0 

5 

0 

Khanchanaburi 

27 

2008-01 


Craseonycteris 

thonglongyai 

- 

- 

7 

0 

6 

0 

Khanchanaburi 

28 

Craseonycteridae cumulation 

- 

10 

0 

11 

0 


2006-03 

Hipposideridae 

Hipposideros 

lekaguli 

1 

0 


- 

2 

0 

Nan 

1 

2006-03 


Hipposideros 

larvatus 

4 

0 


- 

3 

0 

Nan 

1 

2006-03 


Hipposideros 

armiger 

2 

0 


- 

2 

0 

Nan 

1 

2006-04 


Hipposideros 

armiger 

6 

0 

7 

0 

3 

0 

Phayao 

2 

2006-04 


Aselliscus 

stoliczkanus 

1 

0 

1 

0 

1 

0 

Phayao 

2 

2006-04 


Hipposideros 

larvatus 

2 

0 

2 

0 

1 

0 

Phayao 

2 

2006-04 


Hipposideros 

lylei 

2 

0 

2 

0 

2 

0 

Phayao 

2 

2006-04 


Hipposideros 

larvatus 

15 

0 

18 

0 

5 

0 

Chiang Rai 

3 

2006-04 


Aselliscus 

stoliczkanus 

2 

0 

2 

0 

- 


Chiang Rai 

3 

2006-06 


Hipposideros 

armiger 

- 

- 

17 

0 

- 


Ratchaburi 

4 

2006-06 


Hipposideros 

larvatus 

- 

- 

7 

0 

1 

0 

Ratchaburi 

4 

2007-04 


Hipposideros 

armiger 

- 

- 

5 

0 

- 


Uthai Thani 

13 

2006-11 


Hipposideros 

larvatus 

12 

0 


- 

3 

0 

Ratchaburi 

5 

2006-11 


Hipposideros 

armiger 

2 

0 

1 

0 

- 


Ratchaburi 

6 

2006-11 


Hipposideros 

lekaguli 

9 

0 

1 

0 

2 

0 

Ratchaburi 

6 

2007-05 


Hipposideros 

larvatus 

- 

- 

8 

0 

- 


Uthai Thani 

14 

2007-05 


Hipposideros 

pomona 

- 

- 

1 

0 

- 


Uthai Thani 

14 

2007-05 


Hipposideros 

cineraceus 

- 

- 

2 

0 

- 


Uthai Thani 

14 

2007-05 


Hipposideros 

larvatus 

- 

- 

3 

0 

- 


Tak 

15 

2007-05 


Hipposideros 

armiger 

- 

- 

5 

0 

- 


Tak 

16 

2007-05 


Hipposideros 

lylei 

- 

- 

2 

0 

- 


Tak 

16 

2007-09 


Hipposideros 

armiger 

10 

0 

10 


- 


Ratchaburi 

4 

2007-09 


Hipposideros 

larvatus 

1 

0 

1 

0 

- 


Ratchaburi 

4 

2007-03 


Hipposideros 

armiger 

10 

0 

10 

0 

■ 


Chachoengsao 

17 

2007-03 


Hipposideros 

larvatus 

11 

0 



- 


Chachoengsao 

17 

2007-11 


Hipposideros 

pomona 

- 

- 

3 

0 



Nakon Ratchasima 

18 

2007-12 


Hipposideros 

cineraceus 

- 

- 


- 

5 

0 

Khanchanaburi 

24/25 

2007-12 


Hipposideros 

pomona 

- 

- 

2 

0 

1 

0 

Khanchanaburi 

24/25 

2007-12 


Hipposideros 

halophyllus 

- 

- 

2 

0 

4 

0 

Khanchanaburi 

24/25 

2008-09 


Hipposideros 

larvatus 

- 

- 

40 


- 


Chachoengsao 

17 

Hipposideridae cumulation 

90 

0 

163 


35 

0 


2006-03 

Megadermatidae 

Megaderma 

lyra 

2 

0 


- 

1 

0 

Nan 

1 

2006-04 


Megaderma 

spasma 

5 

0 

5 

0 

3 

0 

Lam pang 

12 

Megadermatidae cumulation 

7 

0 

5 

0 

4 

0 


2007-05 

Pteropodidae 

Eonycteris 

spelaea 

- 

- 

2 

0 

- 


Tak 

15 

2007-11 


Rousettus 

leschenaultii 

- 

- 

35 

0 

11 

0 

Lopburi 

21 

2007-11 


Eonycteris 

spelaea 

- 

- 

1 

0 

- 


Nakon Ratchasima 

18 

Pteropodidae cumulation 

- 

38 

0 

11 

0 


2006-03 

Rhinolophidae 

Rhinolophus 

yunanensis 

1 

0 


- 

- 


Nan 

7 

2006-03 


Rhinolophus 

lepidus 

3 

0 


- 

2 

0 

Phitsanulok 

8 

2006-04 


Rhinoloplus 

malayanus 

1 

0 

1 

0 

1 

0 

Phayao 

2 

2006-04 


Rhinoloplus 

sp. 

3 

0 

2 

0 

2 

0 

Phayao 

2 

2007-04 


Rhinolophus 

coelophyllus 

- 

- 

1 

0 

- 


Uthai Thani 

13 

2006-11 


Rhinolophus 

thomasi 

1 

0 


- 

1 

0 

Ratchaburi 

5 

2007-05 


Rhinolophus 

coelophyllus 

- 

- 

6 

0 

- 


Uthai Thani 

14 

2007-05 


Rhinolophus 

malayanus 

- 

- 

1 

0 

- 


Uthai Thani 

14 

2007-05 


Rhinolophus 

thomasi 

- 

- 

1 

0 

- 


Uthai Thani 

14 

2007-05 


Rhinolophus 

coelophyllus 

- 

- 

5 

0 

- 


Tak 

15 

2007-05 


Rhinolophus 

stheno 

- 

- 

1 

0 

- 


Tak 

15 

2007-11 


Rhinolophus 

pearsoni 

- 

- 

5 

0 

- 


Saraburi 

19 

2007-11 


Rhinolophus 

pearsoni 

- 

- 

1 

0 

1 

0 

Lopburi 

20 

2007-12 


Rhinolophus 

malayanus 

- 

- 

1 

0 

- 


Khanchanaburi 

23 

2007-12 


Rhinolophus 

malayanus 

- 

- 

1 

0 

- 


Khanchanaburi 

23 

2007-12 


Rhinolophus 

malayanus 

- 

- 

2 

0 

1 

0 

Khanchanaburi 

24/25 

Rhinolophidae cumulation 

9 

0 

28 

0 

8 

0 


2006-03 

Emballonuridae 

Taphozous 

longimanus 

25 

0 


- 

12 

0 

Lampang 

9 

2006-03 


Taphozous 

longimanus 

12 

0 


- 

4 

0 

Nan 

10 

2006-03 


Taphozous 

longimanus 

9 

0 


- 

1 

0 

Nan 

7 

2006-04 


Taphozous 

longimanus 

22 

0 


- 

- 


Phayao 

2 

2006-04 


Taphozous 

longimanus 

5 

0 

5 

0 

1 

0 

Chiang Rai 

11 

2006-04 


Taphozous 

longimanus 

- 

- 

3 

0 

3 

0 

Phayao 

2 

2007-12 


Taphozous 

longimanus 

2 

0 

2 

0 

- 


Chachoengsao 

22 

2007-03 


Taphozous 

longimanus 

4 

0 

4 

0 

4 

0 

Chachoengsao 

17 

Emballonuridae cumulation 

79 

0 

14 

0 

25 

0 


2006-03 

Vespertilionidae 

Myotis 

siligorensis 

1 

0 


- 

- 


Nan 

7 

2006-04 


Myotis 

siligorensis 

1 

0 


- 

- 


Phayao 

2 

2006-11 


Myotis 

siligorensis 

- 

- 

4 

0 

- 


Ratchaburi 

5 

2006-04 


Myotis 

siligorensis 

1 

0 

1 

0 

- 

0 

Chiang Rai 

3 

2008-01 


Myotis 

siligorensis 

- 

- 


- 

1 

0 

Khanchanaburi 

26 

Vespertilionidae cumulation 

3 

0 

5 

0 

1 

0 


2006-04 

Miniopteridae 

Miniopterus 

schreibersii 

2 

0 

2 

0 

2 

0 

Chiang Rai 

3 


Total number of samples 


All families cumulation 190 


265 


97 


552 

































1694 


M. Ar Couilh et al./Infection, Genetics and Evolution 11 (2011) 1690-1702 


screen another population or species of bats or to assess a putative 
long-lasting infection of a given colony (Table 1). Samples were 
collected from 25 species of bats and were divided and preserved 
in two different media: (i) modified Minimum Essential Medium 
antibiotic supplemented (viral transportation medium) to preserve 
entire virus particles integrity and (ii) RNAlater (RNA preservation 
medium) to limit RNA degradation, both kept on ice before and 
immediately after collection (prior to liquid nitrogen conservation 
in the field). Each sample was immediately aliquoted into four 
tubes (including absolute backup) and kept at -80 °C before RNA 
extraction. Traceability was assured by assigning a unique number 
to each sample. Handling of animals was conducted in accordance 
with the American Society of Mammalogists Guidelines for the 
capture, Handling and Care of Mammals (Gannon and Sikes, 2007). 
Trapping sessions were performed using harp-traps, flip-nets and 
hand-nets. Bats were transferred to individual clean cotton bags 
immediately after being captured and kept in a dark, cool and quiet 
place for one to a few hours. Immediately after sample collection, 
bats were released where they were trapped. Particular attention 
was paid: (i) to avoid disturbance inside the cave and near the col¬ 
ony; (ii) to disinfect material and shoes between sampling sessions 
and study sites; and, (iii) to preserve the surroundings as first 
encountered. 

2.2. RNA extraction 

Qiamp viral RNA mini kit (Qiagen) was used for RNA extraction 
following the manufacturer’s instructions. RNA was eluted in 60 pi 
RNAse-free saline buffer, 7 pi were immediately used for the re¬ 
verse transcription reaction and the remaining 53 pi were ali¬ 
quoted and kept at -80 °C for subsequent detection assays. The 
amount and quality of the extracted RNA was measured using 
Nanodrop and Agilent techniques prior to reverse transcription 
(Mueller et al., 2004; Gallagher and Desjardins, 2006). 

2.3. Reverse transcription (RT), amplification by polymerase chain 
reaction (PCR) and sequencing 

Complement DNA synthesis was performed with M-MLV 
Reverse Transcriptase (Bio-Active) and random hexamer primers. 
Random hexamer primers were used at a final ratio of 0.15 pg 
for 1 pg of RNA. Then, 5 pi of RT products were used in 25 pi final 
volume PCR reaction targeting a 438 nucleotide fragment of the 
RdRp coding region. PCR primers were designed from the RdRp 
coding region by selecting a highly conserved region (nspl2) 
among all Coronavirinae sequences available from online public 
databases in August 2007. Assuming a higher mutation rate for 
the third codon position and according to the published sequence 
variations, degenerate primers were designed: BatCoV pol 15197 
(forward primer: 5'-GGTTGGGAYTAYCCWAARTGTGA-3') and Bat¬ 
CoV pol 15635 (reverse primer: 5'-CCATCRTCMGAHARAATCAT 
CATA-3'), respectively at position 15197 to 15219 and 15635 to 
15614 on Bat SARS-CoV Rfl reference sequence (DQ412042). In or¬ 
der to increase both the sensitivity and the specificity of the detec¬ 
tion, a third primer were designed and combined with the previous 
ones in a semi-nested PCR: BatCoV pol nested 15419 (forward 
nested primer: 5'-GCNAATWSTGTNl Tl AACAT-3'), at the position 
15419 to 15438 on the same reference sequence. The first PCR pro¬ 
duced a fragment of 438 nucleotides and the semi-nested PCR con¬ 
firmed the detection and enhanced significantly the detection limit 
of the semi-nested PCR by amplifying a 216 base pairs fragment 
from the first amplicon. Contaminations were avoided by observ¬ 
ing strict lab procedures, carrying out at each step (mix, PCRI, 
PCR11) in a UV cabinet and using bleached plastic-ware. Manipula¬ 
tions included extraction negative controls at the beginning and at 
the end of each sample series as well as several PCR negative con¬ 


trols per run. The sequencing of the 438 nucleotides PCR products 
was performed in both directions using the external primers. 

2.4. Cell culture and inoculation assay 

Fifty microliters of sample preserved in modified Minimum 
Essential Medium antibiotic supplemented (viral transportation 
medium) were inoculated into 12 wells plates where subconfluent 
cell monolayers were grown. Madin-Darby Canine Kidney (MDCK), 
Vero and Vero clone E6 cells were used in these assays. MDCK were 
chosen because Carnivora are highly suspected to have efficiently 
contributed to the transmission from bats to Humans ( Paguma 
larvata, Nyctereutes procyonides, Melogale moschata ) (Guan et al., 
2003) and several species are susceptible to SARS-CoV infection 
(Martina et al., 2003). The Vero cell line is known to efficiently rep¬ 
licate SARS-CoV (Qinfen et al., 2004). Cells were cultivated in Min¬ 
imum Essential Medium supplemented with 120pg/ml of 
streptomycin, 120 units/ml of penicillin and 10% of Foetal Calf 
Serum (FCS). Once confluence was nearly reached, the supernatant 
was removed and samples were added on the cell monolayer and 
left at 20 °C on a rocking agitator for 60 min. Then, 2 ml of the same 
medium (2% FCS) were added to each well. Plates were then incu¬ 
bated at 37 °C and daily checked on a phase contrast inverted 
microscope for CPE occurrence. 

2.5. Phylogenetic reconstruction 

A 304 or 333 nucleotide region was selected within the se¬ 
quences from both amplicons (Ha AcoV R for 304 or HI Bcov T, 
respectively). The corresponding region from 60 taxa available in 
databases were added to the matrix used in the phylogenetic 
reconstructions. 

1/ Firstly a dynamic homology approach implemented in Poy 
v.4.1.2 was used (Wheeler, 1996; Janies and Wheeler, 2002; Varon 
et al., 2010). This method did not need any prior alignment because 
this step was performed during the cladogram diagnosis and there¬ 
fore established the putative homology statements during the tree 
reconstruction. Poy was used under parsimony as the optimality 
criterion (Wagner algorithm), complemented with the tree bisec¬ 
tion reconnection swapping algorithm (TBR), perturbation injec¬ 
tion, ratcheting and fuse iterations steps (Nixon, 1999). These 
combined steps provided a reasonably fast and accurate way to 
better explore the tree space by escaping from local optima. 
Although RdRp is a conserved gene and the matrix included a great 
number of taxa, resampling taxa tests were used to evaluate 
whether this matrix was subjected to the long branch attraction 
phenomenon which could constitute a limitation of some phyloge¬ 
netic reconstructions. This limitation could be particularly true in 
the case of high evolution rate and very limited number of taxa 
when combined with static homology approach of molecular char¬ 
acters. The over-representation of sub-groups by several similar 
sequences was avoided as it could artificially report changes on 
branches located between internal nodes, leading to the long 
branch attraction issue. 

2/ Secondly, a probabilistic approach (Bayesian) implemented 
in BEAST (Drummond and Rambaut, 2007) was applied to the same 
data set. Sequences were previously aligned using Sea-View (Gouy 
et al., 2010) under the Muscle algorithm (Edgar, 2004) and resulted 
in a frank border matrix showing a unique gap of three codons for 
four sequences (Bulbul, Trush, Munia and Asian leopard CoVs, po¬ 
sition 33 to 36). Then the adequacy of the matrix to 88 substitution 
models was evaluated using jModeltest (Posada, 2008). The three 
codons partition model of evolution was specified (under HYU sub¬ 
stitution model with a gamma site heterogeneity distribution), so 
that each of the three codon position could follow its own evolu¬ 
tion rate. Yule speciation model showed similar estimates as the 
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Table 2 

Asian bat species sampled in coronavirus studies. Some values should be considered cautiously because the related data were not always available and were in some cases 
extrapolated: (*) Details not available, supposedly two samples per bat (respiratory and fecal numbers extrapolated here). ( t ) Expressed as number of bats in original publication, 
details not available. 


Family Genus specie 


PCR Detection 
Respiratory 


Serology Viral ID & Taxonomy (partial) 


2003 / 2004 
Mar 04 

Mar 04 

Mar 04 
July 04 
Dec 04 
Apr 2004 / 
july 2005 
Apr 2004 / 
july 2005 
Nov 04 / 
March 06 
Nov 04 / 
March 06 


HKSAR 

Nanning, Guangxi 

Maoming, 

Guangdong 

Maoming, 

Guangdong 

Nanning, Guangxi 

Nanning, Guangxi 

HKSAR 

HKSAR 

2 Chinese sites 

1 Chinese site 


Hipposid|eridae 


Cynopterus sphinx 
Rousettus leschenaulti 

Rousettus leschenaulti 

Cynoptems sphinx 
Rousettus leschenaulti 
Rousettus leschenaulti 

Cynoptems sphinx 

Rousettus lechenaulti 

Cynoptems sphinx 

Rousettus leschenaulti _ 

Pteropodidae Total: 


42 0 

17 0 


ND 

201 2 1 


Poon et al. 2005 
Li et al. 2005 


Li et al. 2005 
Li et al. 2005 
Li et al. 2005 

Woo et al. 2006 

Woo et al. 2006 

Tang et al. 2006 

Tang et al. 2006 


2003 / 2004 
2003 / 2004 
2003 / 2004 
Apr 2004 / 
july 2005 
Apr 2004 / 
july 2005 
Nov 04 / 
March 06 
Nov 04 / 
March 06 
Nov 04 / 
March 06 
Nov 04 / 
March 06 
Nov 04 / 
March 06 
Nov 04 / 
March 06 


HKSAR 

HKSAR 

HKSAR 

HKSAR 

HKSAR 

13 sites 

1 site 

2 sites 

1 site 

2 sites 
1 site 


Hipposideros armiger 
Hipposideros pomona 
Hipposideros armiger 

Hipposideros armiger 

Hipposideros pomona 

Hipposideros armiger 

Hipposideros larvatus 

Hipposideros pratti 

Hipposideros pomona 

Coelops frithi 

Aselliscus stoliczkanus _ 

Hipposideridae Total: 


Poon et al. 2005 
Poon et al. 2005 
Lau et al. 2005 

Woo et al. 2006 

Woo et al. 2006 

Tang et al. 2006 

Tang et al. 2006 

Tang et al. 2006 

Tang et al. 2006 

Tang et al. 2006 

Tang et al. 2006 


2003/2004 HKSAR 
2003/2004 HKSAR 


2003/2004 HKSAR 
2004 / 2005 HKSAR 


2004/2005 HKSAR 


Miniopteridae 

Miniopterus magnater 
Minioptems pusillus 
Miniopterus schreibersii 
Japanese 

Minioptems magnater 
Minioptems pusillus 


1 


6,25 

26,32 


1 


12,5 

63,16 


Bat-CoV61; Alphacoronavims 
Bat-CoV 61; Alphacoronavims 


3at-CoV 61; Alphacoronavims 


Poon et al. 2005 
Poon et al. 2005 


Poon et al. 2005 
Lau et al. 2005 


Apr 2004 / 
july 2005 

Site B HKSAR 

Minioptems magnater 

102* 


51 

0 

51 


1,96 

ND 


Bat-CoV HKU7b Alphacoronavims 

Woo et al. 2006 

Apr 2004 / 
july 2005 

Apr 2004 / 
july 2005 

Site D HKSAR 

Site E HKSAR 

Minioptems pusillus 

50* 


25 

0 

* 

1 

16 

ND 

ND 


Bat-CoV HKU8-1 Bat-CoV HKU8b 
n=1) Alphacoronavims 

Bat-CoV 61-1/2/3 (as Poon et al. 2005) 
(n=3) Alphacoronavims 

Woo et al. 2006 

Woo et al. 2006 

June 2005 


Minioptems magnater 

25 

■ 0 

ND 


ND 



ND 



Chu et al. 2006 













Bat-CoV 1A-B (8); Bat-CoV HKU8 (1); 

Alphacoronavims & Bat SARS (1) 


Aug 2005 


Minioptems magnater 

44 

m o 

ND 


ND 



ND 


Betacoronavirus 

Chu et al. 2006 

Aug 2005 


Miniopterus pusillus 

1 

0 

ND 


ND 



ND 



Chu et al. 2006 

Dec 2005 

Dec 2005 


Miniopterus magnater 
Miniopterus pusillus 

25 

11 

■ 

ND 

ND 


ND 

ND 



ND 

ND 


Bat-CoV 1A-B Alphacoronavims 

3at-CoV 1A-B Alphacoronavims 

Chu et al. 2006 

Chu et al. 2006 

Mar 2006 


Minioptems magnater 

20 

0 

ND 


ND 



ND 



Chu et al. 2006 

Mar 2006 

2005 

15 sites (3) 

Minioptems pusillus 

Minioptems schreibersi 

10 

270 

0 

ND 

135 

0 

ND 

135 

■ 

12,59 

ND 

ND 


Bat-CoV 1A-B Alphacoronavims 

3tCoV/A773/05 FJ; BtCoV/A911/05 

Chu et al. 2006 

Tang et al. 2006 

Nov 04 

Yichang, Hubei 

Minioptems schreibersi 

1 


ND 


1 

0 


1 

0 


Li et al. 2005 

Miniopteridae Total: 

579 

20 3,45 

297 

7 2,36 

298 

40 

13,42 

1 0 


Vespertilionidae 






Apr 2004 / 
july 2005 

HKSAR 

Myotis chinensis 

16 


8 

0 

8 

0 


ND 



Woo et al. 2006 

Apr 2004 / 
july 2005 

H-HKSAR 

Myotis ricketti 

46 


23 

0 

23 

1 

4,35 

ND 


3at-CoV HKU6 ; Alphacoronavims 

Woo et al. 2006 

Apr 2004 / 
july 2005 

HKSAR 

Nyctalus noctula 

14 


7 

0 

7 

0 


ND 



Woo et al. 2006 

Apr 2004 / 
july 2005 

Sites B-C HKSAR 

Pipistrellus abramus 

28 


14 

■ 

14 

1 

28,57 

ND 


Bat-CoV HKU5-1 site B; Bat-CoV 

HKU5-2,3,5 site C; Betacoronavims 

Woo et al. 2006 

2003 / 2004 

HKSAR 

Pipistrellus abramus 

6 


3 

0 

3 

0 





Poon et al. 2005 

2003 / 2004 

HKSAR 

Myotis myotis 

6 


3 

0 

3 

0 





Poon et al. 2005 

2003 / 2004 

HKSAR 

Myotis ricketti 

10 


5 

0 

5 

0 





Poon et al. 2005 

2004 / 2005 

HKSAR 

Myotis chinensis 

3 


0 


0 



ND 

0 


Lau et al. 2005 

2004 / 2005 

HKSAR 

Myotis hcketti 

2 


0 


0 



ND 

0 


Lau et al. 2005 

2004 / 2005 

HKSAR 

Nyctalus noctula 

2 


0 


0 



ND 

0 


_au et al. 2005 

July 04 

Tianjin 

Myotis ricketti 

42 


21 

0 

21 

0 


ND 



Li et al. 2005 

Nov 04 

Yichang, Hubei 

Myotis altarium 

1 


ND 


1 

0 


1 

0 


Li et al. 2005 

Nov 04 

Yichang, Hubei 

Nyctalus plancyi 

1 


ND 


1 

0 


1 

0 


Li et al. 2005 

2005 

4 sites (1-Hainan) 

Pipistrellus pipistrellus 

54 


27 

0 

27 

■ 

22,22 

ND 


BtCoV/434/05 HA, Betacoronavims 

|Tang et al. 2006 

2005 

8 sites (3-Ahnui, 
Henan, Sichuan) 

Pipstrellus abramus 

82 


41 

0 

41 


34,15 

ND 


BtCoV/355/05 AH, HE, SC, 
Alphacoronavims 

Tang et al. 2006 

2005 

2 sites (1-Hainan) 

Scotophilus kuhlii 

86 


43 

0 

43 

■ 

11,63 

ND 


BtCoV/512/05 HA, Alphacoronavims 

Tang et al. 2006 

Nov 04 / 
March 06 
Nov 04 / 

4 sites 

Myotis daubentonii 

82 


41 

0 

41 

0 


ND 



Tang et al. 2006 

March 06 

1 site 

Myotis mystacinus 

2 


1 

0 

1 

0 


ND 



Tang et al. 2006 

2005 

8 sites (4 Ahnui, 
Yunnan, 
Guangdong, 
Guangxi, Jiangxi) 

Myotis ricketti 

106 


53 

0 

53 

1 

24,53 

ND 


BtCoV/701/05 AH, YN, GD (n=8); 
BtCov/821/05 GX, JX (n=5); 
Alphacoronavims 

Tang et al. 2006 


(continued on next page ) 
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Table 2 ( continued) 


Nov 04 / 
March 06 

2 sites 

Myotis chinensis 

6 

3 

0 

3 

0 


ND 

Nov 04 / 
March 06 
Nov 04 / 

9 sites 

Myotis sp.a 

160 

80 

0 

80 

0 


ND 

March 06 

2 sites 

Nyctalus aviator 

12 

6 

0 

6 

0 


ND 

Nov 04 / 
March 06 

3 sites 

Nyctalus noctula 

34 

17 

0 

17 

0 


ND 

Nov 04 / 
March 06 

1 site 

Scotomanes omatus 

2 

1 

0 

1 

0 


ND 

Nov 04 / 
March 06 

1 site 

Barbastella leucomelas 

2 

1 

0 

1 

0 


ND 

2005 

1 site (1- 
Guangdong) 

Tylonycteris pachypus 

28 

14 

0 

14 

1 

14,29 

ND 

Nov 04 / 
March 06 

1 site 

la io 

16 

8 

0 

8 

0 


ND 

Nov 04 / 
March 06 

1 site 

Murina leucogaster 

10 

5 

0 

5 

0 


ND 

Apr 2004/ 
july 2005 

Sites A, B HKSAR 

Tylonycteris pachypus 

42 

21 

0 

21 

4 

19,05 


Vespertillionidae Total: 

901 

446 

4 0,9 

448 

49 

10,94 


Rhinolophidae 





2003 / 2004 

HKSAR 

Rhinolophus affinis 

4 

2 

0 

2 

0 



2003 / 2004 

HKSAR 

Rhinolophus pusillus 

2 

1 

0 

1 

0 



2003 / 2004 

HKSAR 

Rhinolophus rouxi 

12 

6 

0 

6 

0 



2004 / 2005 

HKSAR 

Rhinolophus affinis 

4 

2 

0 

2 

0 

■ 


ND 0 

2004 / 2005 

HKSAR 

Rhinolophus sinicus 

118 

59 

0 

59 

1 

44,07 

18 1 66,67 

Nov 04 

Yichang, Hubei 

Rhinolophus pusillus 

15 

ND 


15 

0 


ND 

Nov 04 

Yichang, Hubei 

Rhinolophus ferrumequinum 

8 

ND 


8 

■ 

12,5 

4 0 

Nov 04 

Yichang, Hubei 

Rhinolophus macrotis 

11 

3 

0 

8 


12,5 

7 pH 71,43 

Dec 04 

Nanning, Guangxi 

Rhinolophus pearsoni 

41 

11 

0 

30 

■ 

10 

46 H 28,26 

Dec 04 

Nanning, Guangxi 

Rhinolophus pusillus 

8 

2 

0 

6 

0 


6 | 33,33 

Apr 2004/ 
july 2005 


Rhinolophus affinis 

14 

7 

0 

7 

0 


ND 

Apr 2004 / 
july 2005 

Sites D,J HKSAR 

Rhinolophus sinicus 

236 

118 

0 

118 

1 

19,49 

ND 

Apr 2004 / 

Sites F, H, 1 









july 2005 

HKSAR 








ND 

Nov 04 / 
March 06 

27 sites 

Rhinolophus pusillus 

101 

101 

0 

101 

0 


ND 

Nov 04 / 
March 06 

2 sites 

Rhinolophus malayanus 

30 

15 

0 

15 

0 


ND 

Nov 04 / 
March 06 

9 sites 

Rhinolophus affinis 

120 

60 

0 

60 

0 

■ 


ND 


11 sites (3-Ahnui, 
Henan, Sichuan, 









2005 

Hubei, Shandong) 

Rhinolophus ferrumequinum 

82 

41 

0 

41 

■ 

9,76 

ND 


5 sites 

16 sites (1- 

Rhinolophus thomasi 

24 

12 

0 

12 

0 

■ 


ND 

2006 

Shandong) 

10 sites (1- 

Rhinolophus sinicus 

132 

66 

0 

66 

1 

1,52 

ND 

2006 

Shandong) 

Rhinolophus pearsoni 

96 

48 

0 

48 

' 

2,08 

ND 

2004 

11 sites (1-Hubei) 

Rhinolophus macrotis 

76 

38 

0 

38 

1 

2,63 

ND 

Nov 04 / 
March 06 
Nov 04 / 

2 sites 

Rhinolophus rex 

4 

2 

0 

2 

0 


ND 

March 06 

4 sites 

Rhinolophus luctus 

8 

4 

0 

4 

0 


ND 

Nov 04 / 
March 06 

1 site 

Rhinolophus osgoodi 

2 

1 

0 

1 

0 


ND 

Rhinolophidae Total:| 

1148 

599 | 

650 

61 

9,38 



Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 
3tCoV/135/05 GD Betacoronavirus Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 


it-CoV HKU4-1,2,3 Site A; Bat-CoV 

<U4-4, Site B ; Betacoronavirus J: Woo et al. 2006 


Alphacoronavirus 


Poon et al. 2005 
Poon et al. 2005 
Poon et al. 2005 
Lau et al. 2005 


Lau et al. 2005 
Li et al. 2005 
Li et al. 2005 
Li et al. 2005 
Li et al. 2005 
Li ef al. 2005 


Woo et al. 2006 

Bat-CoV HKU2-1, Site D; Bat-CoV 

HKU2-2, Site J (n=2); Alphacoronavirus Woo et al. 2006 


t-SARS-CoV HKU3 (n = 21); 


BtCoV/860/05 SD DQ648838 
BtCoV/970/06 SD DQ648854 
BtCoV/977/06 SD DQ648855 
Alphacoronavirus-, Bt CoV/355 


BtCoV/A1018/06 SD, Betacoronavirus Tang et al. 2006 


| Woo et al. 2006 
Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 


| Tang et al. 2006 
Tang et al. 2006 


BtCoV/970/06 SD ; Alphacoronavirus ; Tang et al. 2006 

7279/04 HB, Betacoronavirus Tang et al. 2006 


Tang et al. 2006 
Tang et al. 2006 
Tang et al. 2006 


constant population size coalescent model. Thus, the Yule specia- 
tion model was chosen as tree prior. The MCMC chain was set to 
10,000,000 states to get an adequate posterior ESS value superior 
to 300. 

The support of the nodes was evaluated via a bootstrap analysis 
(Felsenstein, 1985) for the direct optimization method and poster¬ 
ior probabilities we assessed for the Bayesian method (Fig. 1 and 
SI ). Except for the sequences of Hipposideros larvatus Betacoronavi¬ 
rus Tata (HI BcoV T) and Hipposideros armiger Alphacoronavirus 
Ratcha (Ha AcoV R) which were newly identified in this study, all 
other sequences used in phylogenetic reconstructions were down¬ 
loaded from the National Center for Biotechnology Information 
database (http://www.ncbi.nlm.nih.gov/nucleotide/). Each Gen- 
bank sequence identification number is indicated directly on trees 
(Fig. 1 and SI). 

2.6. Rooting the trees 

In the case of coronavirus phylogeny, rooting with one member 
of the sister group ( Bafinivirus or Torovirus, the two genus of 
Torovirinae sub-family) would have been a good choice to study 
the relation within the in-group Coronavirinae. Unfortunately, the 
region that we used here was too divergent and hardly comparable 
to Torovirinae sequence. However, as Torovirinae sequence can not 


be used to root the trees in this case, we used the results provided 
by a published study (Snijder et al., 2003), that used Torovirinae as 
root. This study, based on an amino acid matrix, inferred a Corona- 
viridae phylogeny, and verified the monophyly of Coronavirus and 
the basal branching of Gammacoronavirus clade. According to these 
conclusions, we chose an avian Coronavirus taxa to root our trees 
and therefore study the Alphacoronavirus and Betacoronavirus gen¬ 
era phylogeny. An additional taxon, from the same clade of the taxa 
used as root ( Gammacoronavirus - Avian Coronavirus Beaudette ), 
was included in the analysis to allow the polarization of the char¬ 
acters (ancestral versus apomorphic) along the tree. 

3. Theory 

3.1. Identifying the phylogenetic signal of filiation in Betacoronavirus 
metagenome 

A genome comprises ribonucleic acid fragments that can inde¬ 
pendently be submitted to various selection pressures of different 
intensity. The history of these fragments can be deeply shaped by 
the function they encode, eventually resulting in incongruent phy¬ 
logenetic signals. The only way to fully recover and clearly exploit 
the phylogenetic signal of a given fragment (by comparison to the 
putative homolog fragment of other taxa) is to analyze it 
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Fig. 2. Colonies of Hipposideros at roost, showing the scattered distribution of individuals, with a zoom on the face of one specimen (bar is equal to 1 cm). (A) Hipposideros 
larvatus (B) Hipposideros armiger. 


independently. Therefore, the combination of several fragments in 
a same analysis should be carefully thought and documented as 
mixing the information they provide can produce a blurred signal 
and therefore can lead to an inaccurate phylogeny. Here we fo¬ 
cused on a fragment of the polymerase sequence only for phyloge¬ 
netic reconstruction of coronaviruses because: (i) the key function 
of this protein involved in vertical transmission and consequently 
in viruses filiation is likely to regulate variation and provide a phy¬ 
logenetic signal that reflects the filiation and, (ii) its relatively short 
length, genome position and associated transcription characteris¬ 
tics make it unlikely to be a putative recombination site. Moreover, 
the analysis of this region by the RDP3 package (Martin et al., 2010) 
did not detect any putative recombination sites. 

During the inference of a species tree, the identification of the 
appropriate vertical transmission of information, is one of the 


classical pitfalls of phylogenetic reconstruction. Therefore, careful 
attention must be paid in the choice of a genetic fragment that re¬ 
flects the species phylogeny and not only the gene history. A given 
genetic fragment, coding for a positively selected function which 
appeared in one species, may benefit another one by horizontal 
transfer (recombination), increasing the fitness of the recombinant 
virus and its progeny. This phenomenon have happened also be¬ 
tween viruses that are actually recognized as members of 
separated families as described for the two independent recombi¬ 
nation events that occurred between a Torovirus ( Coronaviridae , 
Torovirinae ), a Coronavirus and an Orthomyxoviridae (Influenzavirus 
C) (Snijder et al., 1991). These particular events did not involve the 
RdRp sequence. All these reasons lead us to rely on the analysis of 
this RdRp fragment for the phylogenetic reconstruction of 
Coronavirus. 
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3.2. Biodiversity, host phylogeny and micro hot-spots 

According to the distribution of both the known mammalian 
host families, the Rhinolophidae Gray, 1825, and its sister family 
the Hipposideridae Lydekker, 1891, Betacoronavirus- b might be 
widespread over Asia and most of the old world. Moreover, the bio¬ 
diversity of these two mammalian families is surely higher than 
what has been currently described (S.J. Puechmaille, unpublished 
data), which in turn induces a potential greater coronavirus diver¬ 
sity than is currently described. In discordance to what has been 
reported to date and based on the phylogenetic proximity of these 
sister families, Hipposideridae might harbour an equal or richer 
diversity of coronaviruses (at least of the Betacoronavirus genus) 
compared to the Rhinolophidae. Furthermore, the fragmentation 
of the forested habitat, produced by human activities such as land 
use practices and agriculture, leads to the isolation of preserved 
biodiversity spots (surrounded by cultivated fields) that have the 
potential to maintain both host and viral populations, preserving 
by this way the footprints of the natural history of coronaviruses 
in wildlife. 

4. Results 

Twenty-five species of bats belonging to the Rhinolophidae 
(N = 8), Hipposideridae (N = 8), Emballonuridae (N = 1), Molossidae 
(N = 1), Megadermatidae (N = 2), Pteropodidae (N = 2), Miniopteri- 
dae (N = 1) (Table 1 ), Craseonycteridae (N = 1) and Vespertilionidae 
(N= 1) families were sampled. Species were identified using mor¬ 
phological and acoustic characters. 

A total of seven samples among 552 (265 feces, 97 urine collec¬ 
tions and 190 oral swabs), were positive using the RT-PCR on the 
RdRp region described in the material and methods section. This 
amount rose to 28 when the results from the second sampling at 
the study site 17 were included (Table 1 and Fig. IB). All of the po¬ 
sitive samples came from two species: Hipposideros larvatus 
(Fig. 2A) and Hipposideros armiger (Fig. 2B). Two new viruses were 
detected in this study but the isolation trials on cell cultures re¬ 
mained unsuccessful: no cytopathic effect was observed after three 
passages on cell cultures and no viral material was amplified from 
the third passage supernatants with the semi-nested PCR de¬ 
scribed in material and methods. 

4.1. Two coronaviruses from two Hipposideridae: Hipposideros 
armiger Alphacoronavirus Ratcha (Ha AcoV R) and Hipposideros 
larvatus Betacoronavirus Tata (HI BcoV T) 

The Ha ACoV R was detected in bat feces from a unique colony 
(Hipposideros armiger Hodgson, 1835). Among 10 bats tested, four 
were positive (36% - study site 4 - Table 1 ). The BLAST 78% identity 
value indicated that the sequence was related to a bat Alphacorona¬ 
virus reported to infect Rhinolophus ferrumequinum and Rhinolo- 
phus pearsoni in the Chinese Shandong region in 2005 (Table 2) 
and (Tang et al., 2006). For this sequence, 30 bases from the 5' side 
of the PCR products were ambiguous and therefore removed from 
the analysis because the sequencing with one degenerated primer 
gave low quality results. 

Concerning HI BcoV T, the sequences were free of ambiguous 
sites and identical among the molecular isolates amplified from 
feces collected from 22 specimens of bats from the same colony 
(Table 1). The 391 nucleotides long sequences obtained showed 
that this new virus was related to the SARS-CoV clade and shared 
75% nucleotide identity (online BLAST analysis) with the Bat SARSr 
Coronavirus Rft from Rhinolophus ferrumequinum. HI BcoV T was 
detected in fecal samples from one colony of Hipposideros larvatus 
(Hipposideridae Lydekker, 1891), in one sampling site among 26. 


In this colony (estimated to 250 individuals), 3/11 (27%) bats tested 
were positive, then, eighteen months later, 21/40 (52%). This is the 
first evidence of SARSr-CoV long-lasting infection in a bat colony 
(Table 1). 

4.2. Other species and study sites 

Besides the detection of Ha ACoV-R in one site, no other Alpha¬ 
coronavirus was detected from the other study sites. Hipposideros 
armiger individuals positive for Ha ACoV-R shared their roost with 
Hipposideros larvatus individuals, which were negative for the 
detection of Ha ACoV-R. Similarly, besides HI BcoV-T detected in 
one site, no other Betacoronavirus-b was detected from the other 
study sites even if they were populated by the same species as 
the positive site and in two instances (cave 3 and 14 - Table 1), 
with comparable sample sizes. The other species sampled, Taphoz- 
ous longimanus and Hipposideros armiger, that were sharing the 
same roost with Hipposideros larvatus (cave 17 - Table 1), were 
all negative for the detection of HI BcoV T. 

4.3. Phylogenetic reconstruction 

The same global topology was obtained from both direct opti¬ 
mization (Poy) and probabilistic methods (Beast) as apparent dif¬ 
ferences mostly concerned nodes that are not supported by 
respective bootstrap or posterior probability values (inferior to 
95%) (Fig. 1A and SI). Indeed, direct optimization supported Beta¬ 
coronavirus paraphyly but basal nodes exhibited low bootstrap 
support values. Accordingly, the separation between Alphacorona¬ 
virus and Betacoronavirus is not significantly supported in the 
Bayesian phylogeny (posterior probablility value: 0.69 < 0.95). Sev¬ 
eral other nodes remained weakly supported by the data in each 
method. Moreover, the HKU4/HKU5/BtCoV 133 group significantly 
segregated in basal position of the Betacoronavirus in both recon¬ 
structions. Alphacoronavirus exhibited mainly nodes with low sup¬ 
port in both methods of reconstruction despite comparable 
topologies. The Alphacoronavirus detected in this study, Ha AcoV 
R, rooted two alphacoronaviruses reported to infect bats in 
Bulgaria (Drexler et al., 2010). Concerning Betacoronavirus-b topol¬ 
ogy, most of the nodes were well supported. Our RdRp gene based 
phylogenetic reconstructions placed HI BCoV-T as a sister taxon of 
an african Hipposideros Betacoronavirus-b (Fig. 1A and SI). These 
taxa formed together the sister group of all other betacoronavirus- 
es-b, which comprised the monophyletic group of Rhinolophidae 
betacoronaviruses-b. These were rooted by one Hipposideridae 
Betacoronavirus-b: the Australian Bat Coronavirus R. aur, detected 
in Rhinonicteris aurantia. Even paraphyletic, Hipposideridae 
betacoronaviruses mainly grouped together and rooted the Rhinol¬ 
ophidae betacoronaviruses among which is the SARS-CoV (Fig. 1A 
and SI). 

5. Discussion 

The species demarcation criteria proposed to the ICTV by the 
coronavirus study group (http://talk.ictvonline.Org/media/p/ 
1230.aspx) recommend a minimal 90% amino acid sequence iden¬ 
tity threshold over seven domains of the replicase (nsps 3, 5, 12, 
13, 14, 15, 16). Only one domain is available in our data (nspl2) 
but it is located in one of the three most conserved and contiguous 
regions of the replicase (nspll, 12 and 13) (Stadler et al., 2003). 
The relatively low identity values obtained by the BLAST analyses 
(78% for the Ha AcoV-R, and 75% for the HI BcoV-T) in the most 
conserved region of the coronavirus genome suggests that HI BcoV 
T and Ha AcoV R are new species, even though there was only one 
fragment available, instead of the ad hoc set of seven regions pro- 
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posed by the ICTV working group. More investigations are 
necessary to fully characterize these coronaviruses and to get the 
identity values of the seven replicase domains to rule on these 
two taxa description in complete accordance with the ICTV recom¬ 
mendations. For clarity and non-ambiguous taxonomy, these two 
putative new species are proposed to be named Hipposideros larv- 
atus Betacoronavirus Tata (Hi BcoV-T) and Hipposideros armiger 
Alphacoronavirus Ratcha (Ha AcoV R). These Betacoronavirus and 
Alphacoronavirus were detected from the same bat genus Hipposid¬ 
eros (super-family: Rhinolophoidea sensu Teeling et al. (2005), 
family: Hipposideridae Lydekker, 1891). Considering the diversity 
of both Chiroptera and Coronavirinae, we explicitly implied the 
host species and the virus species within the name, thus avoiding 
further possible confusion once other coronaviruses will be discov¬ 
ered and described within these genera (Betacoronavirus and 
Hipposideros ). 

5.1. Sampling pressure and host species involved 

Sampling bias could explain the uneven distribution of corona- 
virus detection among bat families, more specifically concerning 
SARSr-CoV in Rhinolophidae versus Hipposideridae. In contrast 
with other studies, the Hipposideridae family in our study was 
the most sampled family with five times more samples than those 
collected from Rhinolophidae. Only two sampling sites among 28 
(site 17 for HI BcoV-T and site 4 for Ha AcoV-R) explored with var¬ 
ious sampling intensity, gave positive results. For some site, the 
low level of sampling might be insufficient to detect any Coronavir¬ 
inae (Table 1). We cannot exclude that betacoronaviruses are 
hosted by others hipposiderid bat colonies and could circulate 
around Thailand or South-East Asia. Moreover, it is likely that other 
yet unknown betacoronaviruses circulate in Hipposideridae and 
even in other bat families around South-East Asia. An appropriate 
sampling would be able to test these hypotheses. 

5.2. Long-lasting infection in colony, roost-sharing behaviour and host 
phylogenetic affinities are suspected to be species-jumping promoting 
factors 

Chu and collaborators have repeatedly detected alphacoronavi- 
ruses in Miniopterus (Miniopteridae) colonies from Hong Kong. 
Positive detections in bat colonies were obtained over a period 
ranging from 3 months in a Miniopterus pusillus colony to 6 months 
in a Miniopterus magnater colony (Chu et al., 2006). These results, 
suggested that alphacoronaviruses were able to persist through 
time in bat colonies. Another study on Alphacoronavirus (and Beta¬ 
coronavirus - SARSr-Rh-BatCoV -) infected Rhinolophus sinicus bats 
modulated this conclusion. No significant persistence was ob¬ 
served in Rhinolophus sinicus bat colonies and the viral clearance 
was estimated between 2 weeks and 4 months (Lau et al., 2010). 
Moreover, the authors reported a significant loss of body weight 
in bats infected by Betacoronavirus (SARSr-CoVs) that was not ob¬ 
served in Alphacoronavirus infected bats. This loss of weight could 
be the sign of a loss of fitness associated with Betacoronavirus 
infection in Rhinolophus sinicus, without presuming that it would 
be a promoting factor or a consequence of Betacoronavirus infec¬ 
tion. Even though: (i) Betacoronavirus persistence was not observed 
in Rhinolophus sinicus colonies and, (ii) the loss of body mass corre¬ 
lated with Betacoronavirus infection might reduce the fitness of 
Rhinolophus, we report here that a Betacoronavirus was persisting 
for 18 months in a Hipposideridae bat colony. This suggests that 
a Hipposideridae colony might be more tolerant to Betacoronavirus 
than Rhinolophidae over a long period of time or that Betacorona¬ 
virus of Rhinolophidae might have acquired virulence factor limit¬ 
ing long-lasting infection in this host. Further investigations are 
needed to confirm and describe this suspected difference. This 


Betacoronavirus detected so far in Hipposideros larvatus in Thailand 
is restricted to a single cave. Thus, Betacoronavirus long-lasting 
infection in a bat colony might be a rare event. Although rare, 
long-lasting infection can promote cross-species transmission by 
the emergence of new properties in viral populations (Baric et al., 
1999). Despite the coincidence of several promoting factors in 
study site 17 (the long-lasting infection - Table 1 -, the sharing 
of the same small-sized cave by three bat species and the phyloge¬ 
netic proximity of the two Hipposideros species) none of the other 
species were detected positive for Betacoronavirus (Hipposideros 
armiger and Taphozous longimanus were negative - Table 1 ). Inter¬ 
estingly, we noted that the Taphozous longimanus (negative) colony 
was roosting separately from the Hipposideros lan/atus (positive) 
cluster whereas the Hipposideros armiger (negative) colony occu¬ 
pied also another specific area of the cave. Similarly, in site 4 (Ta¬ 
ble 1) the Hipposideros armiger colony, infected by ACoV-Ha, was 
roosting separately from the other (negative) species. Again, this 
typical spatial segregation within a roosting site (e.g. cave) limits 
direct contacts between species and is suspected to account for a 
behavioral barrier to cross-species transmission in the wild, as pre¬ 
viously described for Lentivirus and Felidae (Troyer et al„ 2008; 
VandeWoude et al., 2010). Nevertheless, we mentioned these 
observations here but these mechanisms still remain to be pre¬ 
cisely described in the complex relationship between Rhinolophoi¬ 
dea and betacoronaviruses. 

By contrast, previous studies of Asian coronaviruses that fo¬ 
cused on the species level, suggested that there was no strict asso¬ 
ciation between bat species and Alphacoronavirus species (Chu 
et al., 2006; Tang et al., 2006), and other studies on Betacoronavirus 
in Africa reached similar conclusions (Quan et al., 2010). Species- 
barrier crossing has been reported occasionally for Betacoronavirus 
in Asian, European and African Rhinolophidae or Hipposideridae 
but more frequently in Miniopteridae (Chu et al., 2006; Tang 
et al., 2006; Tong et al., 2009; Drexler et al., 2010). Moreover, Mini¬ 
opterus species exhibit an aggregative behaviour, and mixed spe¬ 
cies colonies are commonly observed as described in literature 
(Benda et al., 2003). Mixed species roosting, inducing direct contact 
between species and therefore promoting transmission, in synergy 
with viral persistence in a colony and the phylogenetic proximity 
of the hosts, is likely to promote cross species transmission of 
alphacoronaviruses. As a result, a non-specific host-virus associa¬ 
tion for Alphacoronavirus in Miniopterus follows the same model 
as what has been described for these bats and their ectoparasites 
(Bruyndonckx et al., 2009). The apparent higher frequency of 
cross-species transmission (or the absence of strict virus-host asso¬ 
ciation) in Miniopteridae compared to Hipposideridae and Rhinol- 
phidae might partially be explained by their different roosting 
behaviour. 

5.3. Phytogeny and biogeography: Relation between betacoronaviruses 
filiation and host distribution 

The two methods of phylogenetic reconstructions employed 
here supported comparable relationships within Coronavirinae 
(Fig. 1 and SI ). The uncertainty associated with some relationships 
was often similarly reported by each method. These phylogenetic 
issues might be reduced by increasing the size of the analysed frag¬ 
ment, that would increase the number of informative characters 
and therefore provide more information to support nodes. Adding 
new taxa would also improve the phylogeny by providing probable 
currently lacking information, within or around groups that are not 
yet well supported: (i) Betacoronavirus sub-groups relationships 
might be improved by looking for new taxa related to the group 
HKU9/Bat CoV Philippine Dilim (weak posterior probability: 0.59 
and bootstrap value: 0.24 are reported for this node); (ii) Other se¬ 
quences from new Betacoronavirus (or taxa from a new group) 
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might help to clarify the relationship between Betacoronavirus, 
Alphacoronavirus and putative new groups as discordance is ob¬ 
served between phylogenies and none of the method provide sig¬ 
nificant support to the nodes that link Betacoronavirus and 
Alphacoronavirus (posterior probability: 0.69 and bootstrap values: 
0.21; 0.06; 0.24). 

Interestingly, our study revealed that the Betacoronavirus de¬ 
tected in Hipposideridae in South-East Asia was related to those 
that currently circulate in Africa, according to the analysis of a 
RdRp fragment. Our analyses strengthen the diversity of Betacoro- 
navirus- b clade and extend their broad distribution. Indeed, SARS- 
CoV emerged from this Betacoronavirus- b clade (Fig. 1A) which is 
present in Africa, Europe, Asia and even Australia. 

Quan and collaborators pushed forward the possible occurrence 
of migration events from Africa to Asia to explain the basal phylo¬ 
genetic position of an African Betacoronavirus (compared to Asian 
rhinolophidae SARSr-CoV) together with the apparent absence of 
SARSr-CoV in Asian Hipposideridae (Quan, 2010). With the detec¬ 
tion of Betacoronavirus in hipposiderids in Asia, this hypothesis, 
while still possible, becomes less necessary to explain phylogenetic 
and geographic patterns of Hipposideridae and Rhinolophidae 
betacoronaviruses. The clustering of Betacoronavirus- b that corre¬ 
lates with the host phylogeny rather than with virus geolocaliza¬ 
tion suggests that the phylogenetic relationships between 
betacoronaviruses-b is mainly driven by the host phylogeny 
(Fig. 1 and SI). Accordingly, the phylogenetic relationships of the 
HI BCoV detected here (from South-East Asia) with the SARSr- 
CoV (from Africa) and other betacoronaviruses support the follow¬ 
ing hypothesis. Betacoronaviruses-b ancestors, meaning 
SARSr-CoVs ancestors, could have been historically hosted by the 
common ancestor of the Rhinolophidae and Hipposideridae and 
could have later evolved independently in the lineages leading to¬ 
wards Rhinolophidae and Hipposideridae betacoronaviruses. This 
point of view is supported by the topology of Betacoronavirus clade 
that exhibits the monophyly of Rhinolophidae betacoronaviruses 
despite their broad geographic distribution (from Africa and cen¬ 
tral Europe to South-East Asia and Northern Australia - Fig. 1) 
and shows the close relationship of Hipposideridae betacoronavi¬ 
ruses. Nevertheless, it is unlikely that putative host switching be¬ 
tween Rhinolophidae and Hipposideridae betacoronaviruses 
(ancestors), eventually associated with ancient migration events, 
would still be detectable. More investigations are necessary to ad¬ 
dress these points and test corresponding hypotheses. As well, new 
investigations would be necessary to confirm this suggested pre¬ 
dominant influence of the host phylogeny on Betacoronavirus- b 
clustering and whether other lineage would exist in Hipposideri¬ 
dae, as suggested by the report of a Betacoronavirus from Australia 
(C. Smith pers. com. - Fig. 1 ). 

The position of the Hipposideridae Betacoronavirus from 
Australia (basal to the Rhinolophidae betacoronaviruses), suggests 
the existence of an intermediary subgroup (Fig. 1 and SI) between 
Hipposideridae and Rhinolophidae betacoronaviruses that sur¬ 
vived within a rare endemic Australian bat species. The phyloge¬ 
netic status of this host, Rhinonicteris aurantia, is not well defined 
but this species is described as intermediary between Hipposider¬ 
idae and Rhinolophidae, while considered an Hipposideridae 
(Hand and Kirsch, 1998; Simmons, 2007; Armstrong, 2006; Arm¬ 
strong com. pers.). Albeit this putative intermediary betacoronavi¬ 
rus group count a unique member to date and its associated 
posterior probability is weak, its phylogenetic position correlates 
with that of its host. Therefore, we are looking forward new inves¬ 
tigations that would address these elements and screen betacoro¬ 
naviruses in other Hipposideridae and Rhinolophidae. 

The ecological niche of these bats and their behaviour provide 
them with frequent opportunities to spread viral variants among 
other nocturnal mammals such as other Chiroptera or Carnivora 


Viverridae (civet family) (Song et al., 2005; Guan et al., 2003; 
Gonzalez et al., 2008). Many factors have been suggested to explain 
the occurrence of emerging viruses in bats such as diversity (in 
term of species and ecology), ability to fly, immunology, physiol¬ 
ogy, seasonal migration and roosting behaviour (Dobson, 2005; 
Kunz and Fenton, 2005; Calisher et al., 2006). Another stated that 
population genetic structure correlates with viral richness (Turm- 
elle and Olival, 2009). In the case of Betacoronavirus, the role of: 
(i) the fragmented biotope observed today (the Thai forest biotope 
has drastically reduced over the past decades, leading to the forma¬ 
tion of species refuge areas promoting species interactions (Harvey 
et al„ 2006; Trisurat et al., 2010)) and, (ii) the segregation behav¬ 
iour and the mobility patterns of bats (metapopulations, migra¬ 
tions), should be investigated for their impact on the mosaic 
pattern of Betacoronavirus. 

6. Conclusion 

This study suggests that Hipposideridae have been underesti¬ 
mated in the study of coronaviruses in the wild. Both an Alphacor¬ 
onavirus and a Betacoronavirus (more precisely the sub-group b 
that includes SARS-CoV) were detected in Hipposideridae in 
South-East Asia (Thailand). In contrast to what is reported for Rhi¬ 
nolophidae, Hipposideridae are able to host Betacoronavirus close 
to SARS-CoV over long periods of time, suggesting that they could 
be a source of virus for other species under specific circumstances. 
Moreover, the phylogeny of the RdRp, when compared to the 
biogeography, suggests that these two bat families may share an 
ancient relationship with independent lineage of Betacoronavirus. 
Consequently, and by contrast to what have been generally re¬ 
ported for coronaviruses, we suggest that the host phylogeny glob¬ 
ally drive the pathogen clustering in the particular relationship 
between Betacoronavirus and Rhinolophoidea. As for other emerg¬ 
ing viruses, interactions between wildlife refuges (whether consid¬ 
ered in spatial or organic dimensions) and human populations, 
likely enhanced by fragmentation, changes, cultural behaviour 
and secular traditions, contribute to drive the species barrier cross¬ 
ing of coronaviruses to humans. 
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